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Abstract 

We determine parameters for MSW and vacuum oscillations (active and ster- 
ile neutrinos) that are allowed by separate, and collective, imposition of the 
constraints from total event rates in the chlorine, GALLEX, SAGE, and Su- 
perKamiokande experiments (504 days), the SuperKamiokande energy spec- 
trum, and the SuperKamiokande zenith-angle dependence. The small mixing 
angle MSW solution is acceptable at 7% C.L. (8% for sterile u's) and the 
vacuum solution is acceptable at 6% C.L. . The best-fit global MSW solution 
for active neutrinos is: Am 2 = 5 x 10~ 6 eV 2 , sin 2 29 = 5.5 x 10~ 3 (and for 
sterile neutrinos: Am 2 = 4 x 10~ 6 eV 2 , sin 2 29 = 7 x 10~ 3 ). For vacuum os- 
cillations, the best-fit solution is: Am 2 = 6.5 x lCr n eV 2 , sin 2 26 = 0.75 . An 
arbitrary combination of undistorted (no oscillations) pp, 7 Be, 8 B, and CNO 
neutrino fluxes is inconsistent with the combined data sets at the 3. 5a C.L., 
independent of astrophysical considerations. We use improved calculations of 
solar model fluxes, neutrino absorption cross sections and energy spectra, and 
a detailed evaluation of regeneration effects. 
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I. INTRODUCTION 



It is now 30 years since the first recognition of the solar neutrino problem K-H]- In the 
first two decades of solar neutrino research [[| , the problem consisted only of the discrepancy 
between theoretical calculations based upon on a standard solar model (with the implicit 
assumption that neutrinos created in the solar interior reach the earth unchanged) and the 
observations of the capture rate in the chlorine solar neutrino experiment. 

In recent years, four new experiments (Kamiokande [EJ, GALLEX 0, SAGE and 
Super Kamiokande [^-[T0[|) have reported results. All four experiments confirm the original 
detection of solar neutrinos with lower neutrino fluxes than predicted by standard solar mod- 
els. In addition, the Kamiokande and SuperKamiokande experiments demonstrate directly 
that the neutrinos come from the sun by showing that recoil electrons are scattered in the 
direction along the sun-earth axis. 

In April 1996, the SuperKamiokande experiment initiated a new era of high-precision and 
high-statistics solar neutrino research. The first 504 days of data from SuperKamiokande 
|T0| , when combined with data from earlier experiments on solar neutrinos, provide impor- 
tant constraints on the MSW HIT! and vacuum oscillation fOl solutions of the solar neutrino 



problem. 

The definitive analysis of the implications of the SuperKamiokande data set must await 
the comprehensive Monte Carlo study that can only be performed by the SuperKamiokande 
Collaboration. However, the excellent agreement of the results from the first 300 days of 
SuperKamiokande operation || with the results obtained after 374 days || and after 504 
days fTJJ shows the robustness of the results. Therefore, with more than a year's worth 
of data available, this is an appropriate time to take stock of what has been achieved and 
what further challenges lie ahead. The analysis presented here will, we hope, be useful in 
guiding experimental plans for other detectors such as SNO |13[ and BOREXINO WM. We 



also hope that our summaries of the analysis techniques and the theoretical input data (see 
especially Sec. |TJ and the Appendix), as well as the indicated results, will be helpful to 
others who will make similar studies. The results from different theoretical analyses should 
be compared with each other, and with the comprehensive studies by the SuperKamiokande 
Collaboration using their detailed Monte Carlo simulation, in order to test the robustness 
of the inferences about neutrino parameters. 

In this paper, we explore the implications of the 504 day data set from SuperKamiokande 
together with the results from the chlorine, GALLEX, and SAGE experiments. For a concise 
summary of our conclusions, the reader is advised to skip directly to Sec. |V111| and then to 
return to this introduction. 

We use improved neutrino interaction cross sections JT5]HTE|] , the most accurate neutrino 
spectra 0,0, and the results of a recent reevaluation, the BP98 model, by Bahcall and 
Pinsonneault of the standard model neutrino fluxes [2C|. The BP98 model is based upon a 
comprehensive examination of all of the available nuclear fusion data |21] that was carried out 
under the auspices of the Institute of Nuclear Theory (INT). For solar neutrino research, the 
most important nuclear physics parameter is the low energy cross section factor, Sn, for the 
reaction 7 Be(p, 7) 8 B, which gives rise to the critical 8 B neutrinos. The INT normalization, 
MINT) 

SMGLT) 



19-2 e V b [BTJ, is about la less than the previously standard Caltech (CIT) cross 



section factor [22 



(22.4 ± 2.1) eV b. The CIT normalization was computed 
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by taking the weighted average over all the published experimental data while the INT 
normalization was computed by including only the results from the two most recent and 
best documented experiments. Because the uncertainty in the 8 B production cross section 
is the most important uncertainty in predicting solar neutrino fluxes ||20|| , we present results 
in this paper for both the INT and the CIT normalizations of Su. 

We also present calculations for oscillations into sterile neutrinos as well as into the more 
familiar active neutrinos. In what follows, we shall always mean oscillations into active 
neutrinos unless we explicitly include the adjective 'sterile.' We do not consider here the 
intermediate case of oscillations partly into active neutrinos and partly into sterile neutrinos, 



although this is a logical possibility. For plausible assumptions (see e.g. |23[), consistent 
with the results of the CHOOZ experiment ||24|| , the oscillations of solar neutrinos are well 
described by the two-generation formalism. Therefore, we consider explicitly here only two 
neutrino generations. 

We begin by considering in Sec. [TT] how well (or rather how poorly) the results of solar 
neutrino experiments are described by the combined predictions of the standard solar model 
and the minimal electroweak theory (which implies that nothing happens to the neutrinos 
after they are created). In Sec. |T| we briefly summarize the ingredients and techniques used 
in our analysis and give references to the original sources for the improved neutrino flux cal- 
culations, the associated uncertainties in the fluxes, the improved neutrino energy spectra 
and interaction cross sections, the methods including theoretical errors, and the techniques 
for carrying out the theoretical calculations. Details of the statistical analysis are provided 
in the Appendix. We determine in Sec. |IV| the regions that are allowed in neutrino param- 
eter space for MSW and vacuum oscillations to either active or sterile neutrinos provided 
only the total event rates in the neutrino experiments are considered. In Sec. we deter- 
mine the implications of the zenith-angle dependence of the SuperKamiokande event rates. 
We quantify in Sec. [VI| the distortion of the recoil electron energy spectrum measured by 
SuperKamiokande, determining the slope parameter and the excluded regions of oscillation 
parameters. We impose in Sec. [VIl] all of the constraints, total rates, electron recoil energy 
spectrum, and zenith-angle dependence, in global fits and determine the range of oscillation 
parameters that are consistent with all the data. We present in Sec. |V1U| our summary and 
overview of where we stand in understand the discrepancies between standard model pre- 
dictions and the results of solar neutrino experiments and in the determination of neutrino 
parameters from solar neutrino experiments. 



II. WHAT IS ALL THE FUSS ABOUT? 

Why are so many papers being written about non-standard physics implied by solar 



neutrino experiments? This section provides two answers to this question. In Sec. [II A[ we 
show that all of the 19 standard solar model calculations published in refereed journals in 
the last 10 years predict neutrino fluxes that are in reasonable agreement with each other. 
In Sec. irBl, we show that the measured rates and their uncertainties in solar neutrino 



experiments are inconsistent with any combination of the solar neutrino fluxes that does 
not include a spectrum distortion-which requires physics beyond the standard electroweak 
model. We summarize in Sec. [II C| the comparisons between the predictions of the stan- 
dard model-minimal electroweak theory plus standard solar model-and the results of solar 
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neutrino experiments. 



A. The last decade of standard solar models 



Figure [I] displays the calculated 7 Be and 8 B neutrino fluxes for all 19 standard solar 
models with which we are familiar which have been published in the last 10 years in refereed 
science journals. The fluxes are normalized by dividing each published value by the flux 
from the BP98 solar model f20fl ; the abscissa is the normalized 8 B flux and the ordinate is 



the normalized 7 Be neutrino flux. The rectangular box shows the estimated 3a uncertainties 
in the predictions of the BP98 solar model. The abbreviations, which indicate references to 
individual models, are identified in the caption of Figure [l]. 

All of the solar model results from different groups fall within the estimated 3a uncer- 
tainties in the model predictions (with the exception of the Dar-Shaviv model whose results 
have not been reproduced by other groups). This agreement demonstrates the robustness of 
the predictions since the calculations use different computer codes (which achieve varying 
degrees of precision) and involve a variety of choices for the nuclear parameters, the equation 
of state, the stellar radiative opacity, the initial heavy element abundances, and the physical 
processes that are included. 

The largest contributions to the dispersion in values in Figure [I] are due to the choice 
of the normalization for S17 (the production cross-section factor for 8 B neutrinos) and the 
inclusion, or non-inclusion, of element diffusion in the stellar evolution codes. The effect 
in the plane of Fig. [I] of the normalization of Sn is shown by the difference between the 
point for BP98 (1.0,1.0), which was computed using the INT normalization, and the point 
at (1.18,1.0) which corresponds to the BP98 result with the CIT normalization. 

Helioseismological observations have shown recently |25| that diffusion is occurring and 



must be included in solar models, so that the most recent models shown in Fig. |l] now all 
include helium and heavy element diffusion. By comparing a large number of earlier models, 
it was shown that all published standard solar models give the same results for solar neutrino 
fluxes to an accuracy of better than 10% if the same input parameters and physical processes 
are included [^,0 • 

How do the observations from the solar neutrino experiments agree with the solar model 
calculation? 

Table | summarizes the solar neutrino experimental rates that have been measured in the 
five experiments. We have compared the observed rates with the calculated, standard model 
values, combining quadratically the theoretical solar model and experimental uncertainties, 
as well as the uncertainties in the neutrino cross sections. Since the GALLEX and SAGE 
experiments measure the same quantity, we treat the weighted average rate in gallium as one 
experimental number. We adopt the SuperKamiokande measurement as the most precise 
direct determination of the higher-energy 8 B neutrino flux. 

Using the predicted fluxes from the BP98 model, the \ 2 f° r the fit to the three experi- 
mental rates (chlorine, gallium, and SuperKamiokande) is 

Xssm(3 experimental rates) = 61 . (1) 

The result given in Eq. (|I|), which is approximately equivalent to a 20a discrepancy, is 
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a quantitative expression of the fact that the standard model predictions do not fit the 
observed solar neutrino measurements. 



B. Model-independent tests 



Suppose (following the precepts of Hata et al. [^8) , Parke |2{J , and Heeger and Robert 



son 



|5D| ]) we now ignore everything we have learned about solar models over the last 35 years 
and allow the important pp, 7 Be, and 8 B fluxes to take on any non-negative values. What 
is the minimum value of x 2 for the 3 experiments, when the only constraint on the fluxes 
is the requirement that the luminosity of the sun be supplied by nuclear fusion reactions 
among light elements? We include the nuclear physics inequalities between neutrino fluxes 



(see section 4 of Ref. |JT|) that are associated with the luminosity constraint and maintain 
the standard value for the essentially model-independent ratio of pep to pp neutrino fluxes. 



1. With SSM value for CNO neutrinos 

We begin by allowing the pp, 7 Be, and 8 B neutrino fluxes to be arbitrary parameters 
(subject to the luminosity constraint), but constrain the 13 N and 15 fluxes to be equal 
to the values predicted by the standard solar model [pOR . There is therefore one degree of 



freedom for three experiments plus the luminosity constraint and three freely chosen neutrino 
fluxes. The best fit for arbitrary pp, 7 Be, and 8 B fluxes is obtained for 7 Be/( 7 Be) SSM = 0.0, 
8 B/( 8 B) SSM = 0.44, and pp = 1.08, where 

Xminimum(3 experimental rates; arbitrary pp, 7 Be, 8 B) = 24.1 . (2) 

The best-fit solution differ considerably from the measured values for the radiochemical 
experiments, with 3.4 SNU for chlorine and 95.5 SNU for gallium, but is in good agreement 
with the measured value (0.474 of BP98) for SuperKamiokande (cf. Table 0). 
There are no acceptable fits at a C.L. of more than 99.99% (4a result). 



2. CNO neutrinos assumed missing 

The fit can be improved if we set the CNO neutrino fluxes equal to zero. f\ Then, the same 
search for arbitrary pp, 7 Be, and 8 B neutrino fluxes leads to a best fit with 7 Be/( 7 Be) SSM = 
0.0, 8 B/( 8 B) sgM = 0.46, and pp/(pp)ssM = 1-10 .. The minimum value of \ 2 is 

Xmmimum(3 experimental data; arbitrary pp, 7 Be, 8 B; CNO = 0) = 7.3 . (3) 

The best-fit solution has 2.9 SNU for chlorine and 85.3 SNU for gallium, and 0.46 of BP98 
for SuperKamiokande (cf. Table |). Although these three best-fit values are not far from the 



*If the CNO neutrino fluxes are allowed to vary as free parameters, the minimum x 2 is achieved 
for zero CNO fluxes. 
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three measured values, the best-fit values were found by an extensive computer search with 
three free parameters and only one physical constraint, that the nuclear energy generation 
rate correspond to the total luminosity of the sun. 

There are no acceptable solutions at the 99% C. L. (~ 3a result). 

Figure [I] shows the best-fit solution and the lcr -3a contours. The la and 3a limits were 
obtained by requiring that \ 2 — Xmin + ^X 2 ; where for la 5\ 2 = 1 and for 3a S\ 2 = 9. All 
of the standard model solutions lie far from the best-fit solution and even lie far from the 
3a contour. 

C. Discussion of comparisons with standard model 

The searches for best-fit solutions that are described by Eq. (0) and Eq. (|) have three 
independent experimental data points and 1 additional constraint (the luminosity constraint) 
with three free parameters (the pp, 7 Be, and 8 B neutrino fluxes). We conclude from Eq. (Q) 
and Eq. @ that all solutions with undistorted energy spectra are ruled out at the 99% C.L. 
(with the CNO neutrino fluxes set equal to zero, see Eq. ([$]) or at more than 99.9% C.L. 
(with the CNO fluxes equal to the values predicted by the standard solar model). 

The results in this section (see Eq. |T] and Fig. [I]) show that the combined standard model, 
the standard solar model plus minimal standard electroweak theory, provides a bad fit to the 
observed rates in solar neutrino experiments. Moreover, Eq. ([|) and Eq. @ show that one 
cannot get a good fit to the observed rates using any combination of undistorted neutrino 
energy spectra. All so-called astrophysical solutions give a poor description of the observed 
experimental rates. This is what all the fuss is about. 

If we drop the physical requirement that the fluxes be positive definite, the minimum 
X 2 (Xmin = 0-1) occurs-quite remarkably- for a negative value of the 7 Be flux; 0( 7 Be) m i n = 
—0.45 x 0( 7 Be)ssM- This unphysical result is a reflection of what has become known as the 
problem of " the missing 7 Be solar neutrinos.". One reason that the 7 Be neutrinos appear 
to be missing (or have a negative flux) is that the two gallium experiments, GALLEX and 
SAGE, have an average event rate of 72.3 ± 5.6 SNU, which is fully accounted for in the 
standard solar model by the fundamental pp and pep neutrinos (72.4 SNU) p6| , p7| , p0| . In 
addition, the 8 B neutrinos that are observed in the Kamiokande and SuperKamiokande 
experiments will produce about 6 SNU in the gallium experiments, unless new particle 
physics changes the neutrino energy spectrum. A second reason that the 7 Be flux appears 
to be missing is that (see pUfl ) the SuperKamiokande 8 B neutrino flux alone corresponds 
to 2.70(1 ± 0.11) SNU in a chlorine detector (combining quadratically the theoretical cross 
section errors with the SuperKamiokande measurement errors), which is to be compared to 
the total observed rate in the Homestake chlorine experiment of (2.56 ± 0.23) SNU. The 
observed rate includes neutrinos from 8 B, 7 Be, CNO, and pep. Obviously, there is no room 
in the measured chlorine rate for a significant 7 Be neutrino contribution (expected to be 
1.15 SNU based upon the standard solar model p0f). 

For all these reasons, we will consider in the remainder of this paper theories in which 
neutrino oscillations change the shape of the neutrino energy spectra. 
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III. INGREDIENTS AND TECHNIQUES 



Many authors have reported the results of refined studies of the MSW P^|-^T[] and the 
vacuum |57],^,[1(],[|^] solutions of the solar neutrino problems. The techniques for this analy- 
sis are therefore well documented in the literature and we only note briefly here those aspects 
of the calculation that are often not treated in the optimal manner in the literature or for 
which less accurate data are sometimes used. Our \ 2 analysis of the data follows closely 
the prescriptions in P3|j3Tf . We adopt the procedures of Fogli and Lisi (see Ref. J33|) in 
including theoretical errors. The uncertainties in the input model parameters that influence 

2D| . f27| ]. We use the improved neutrino interac- 



the neutrino fluxes are taken from Refs. 
tion cross sections for each detector given in Refs 



15 18 1 and the neutrino spectra given 



in Refs. 



We include the published energy resolution and trigger efficiency of the 



Kamiokande detector || and SuperKamiokande |8|-[T0|. For the MSW solution, we use the 
analytical description of the neutrino survival probabilities from [33] which allows the aver- 
aging over the neutrino production regions and the neutrino spectra to be done accurately 
with a reasonable amount of computer time. 

We obtain allowed regions in Am 2 - sin 2 29 parameter space by finding the minimum \ 2 
and plotting contours of constant x 2 — Xmin + <^X 2 where A% 2 = 5.99 for 95% C.L. and 
9.21 for 99% C.L. . In all the figures in this paper, we show results at the 99% C.L. 
We describe the statistical analysis in detail in the Appendix. 



IV. FITS TO THE AVERAGE EVENT RATES 

In this section, we determine the allowed range of oscillation solutions using only the total 
event rates in the ClAr (Homestake), GALLEX, SAGE, and SuperKamiokande experiments.^ 

The average event rates in these four solar neutrino experiments, summarized in Table |, 
are robust and seem unlikely to change significantly. The results of the chlorine experiment 
have been summarized in detail recently ||45|| ; the measured value has been stable for two 
decades and is known relatively precisely. The Kamiokande and SuperKamiokande [@-[l(J 
experiments are in agreement to within la and the precision of the SuperKamiokande ex- 
periment is now very high. We have used the more precise value of the SuperKamiokande 
experiment to represent the total rate above 6.5 MeV in the water Cherenkov experiments. 
Both the GALLEX || and the SAGE [0] experiments have yielded accurate measurements 
and the efficiency of both detectors has been tested with 51 Cr sources. We have used the 
weighted average of the measured rates in the two gallium detectors. 



We present in Sec. [IV A| the allowed range of solutions involving MSW oscillation into 
active or sterile neutrinos. In Sec. [IV B| , we present the corresponding results for vacuum 
oscillations. Sec. |IV C| describes the dependence of the inferred oscillation parameters on 



^We have carried out identical calculations including both the Kamiokande and the Su- 
perKamiokande rates. The results are essentially unchanged from what we find including only 
the SuperKamiokande rate, since the quoted uncertainty in the Kamiokande rate is much larger 
than the uncertainty in the SuperKamiokande rate. 
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the most uncertain input parameter, the low-energy cross section factor for the 7 Be(p, 7) 8 B 
reaction. We show in Sec. [VP that energy-independent oscillations are unacceptable at the 
99.8% CL. . Finally, we summarize in Sec. |1V E| our results on the analysis of the average 
event rates. 

The calculation of the allowed range of predicted 7 Be flux (to be measured in the BOREX- 
INO experiment flip, the demonstration that a constant suppression factor is a disfavored 
description, and the evaluation of the dependence of the inferred neutrino parameters on 
the low-energy cross section factor for the 7 Be(p, 7) 8 B reaction are special features of this 
section. 



A. MSW solutions 

1. Active neutrinos 

The best fit is obtained for the small mixing angle (SMA) solution: 

Am 2 = 5.4 x 10" 6 eV 2 , (4a) 
sin 2 20 = 6.0 x 10" 3 , (4b) 

which has Xmm = 1-7- There are two more local minima of x 2 ■ The best fit for the well 
known large mixing angle (LMA) solution occurs at 

Am 2 = 1.8 x l(T 5 eV 2 , (5a) 
sin 2 29 = 0.76, (5b) 

with Xmin = 4-3. There is also a less probable solution, which we refer to as the LOW 



solution (low probability, low mass), at 

Am 2 = 7.9 x 10" 8 eV 2 , (6a) 
sin 2 29 = 0.96. (6b) 

with Xmm = 7-3. The LOW solution is acceptable at the 99% C.L., but is not acceptable at 
the 95% C. L. . To find an appreciable probability for the LOW solution, one must include 
the regeneration effect |43[ . 



How do the results given in Eq. (|)-Eq. (|(J) differ from our pre-SuperKamiokande (1997) 
analysis? (The 1997 study, Ref . |43| , was carried out before either the Bahcall-Pinsonneault 



98 solar model or the SuperKamiokande experimental results were available.) There are no 
very large changes in the best-fit values for either the mass differences or the mixing angles 
of the three solutions, although the x 2 fits are less good (but well within the statistical 
uncertainties, cf. Eqs. (l)-(3) of Ref. However, the best-fit neutrino parameters are 

shifted towards smaller mixing angles for the SMA and towards larger mixing angles for the 
LMA and LOW solutions. For the SMA solution, this has the important result that the 
expected Day-Night asymmetry is greatly decreased (see Sec. |V|) 

Figure ^| shows the 99% CL. allowed regions in the plane defined by Am 2 and sin 2 29. 
The black dots within each allowed region indicate the position of the local best-fit point 
in parameter space. The results shown in Fig. [| were calculated using the predictions of 
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the 1998 standard solar model of Bahcall and Pinsonneault |2(J , which includes helium and 



heavy element diffusion and uses the recent reevaluation of solar fusion cross sections ||21|| ; 
the shape of the allowed contours depends only slightly upon the assumed solar model (see 
Fig. 1 of Ref. 0). 

The BOREXINO experiment will measure the v — e scattering rate for the 0.86 MeV 7 Be 
line. We have calculated the allowed range for the scattering rate at 99% CL. for the different 
MSW solutions discussed above. The rates are given in the following equation relative to 
the 1998 Bahcall and Pinsonneault standard model [EOfl are: 



(0Cr) 7 Be SMA _ n 9 o+0.24 / 7q \ 
{(Per) BP98 



= o.59l°;- , (7b) 



^Tge LOW _ n c-o+0.06 



BP98 



0.59^ . (7c) 



2. Nuclear physics uncertainties 

After more than 35 years of progressively more accurate measurements of input param- 
eters and more precise solar modeling, the largest recognized uncertainties that afflict the 
prediction of solar neutrino fluxes are associated with the laboratory cross sections at low 
energies for the crucial 3 He(a,7) 7 Be and the 7 Be(p, 7) 8 B reactions ||20|| . The uncertainty 



in the 3 He — 4 He cross section is 9.4%, la ||21|| . The uncertainty in the 7 Be + p reaction is 
asymmetric and is, on average, 10.6% ||21||. The flux of 8 B neutrinos is directly proportional 



to the rate of the 7 Be + p reaction and the flux of the 7 Be neutrinos is approximately linearly 
proportional to the rate of the 3 He — 4 He reaction. 

How much do the uncertainties in the nuclear physics parameters affect the accuracy 
with which one can determine neutrino parameters? In order to answer this question, 
we have calculated the allowed regions for MSW solutions, analogous to those considered 
in Sec. |1V A 1| , but we assumed in the present case either that the uncertainty in the 



3 He(a,7) 7 Be or the 7 Be(p, 7) 8 B reaction was equal to zero. 

Figure | shows that the MSW allowed regions are reduced only slightly when the uncer- 
tainties associated with either the 3 He(a,7) 7 Be (cross-section factor S34) or the 7 Be(p, j) 8 B 
reaction (cross-section factor Sn) are artificially decreased to zero. Comparing Fig. |3] and 
Fig. 0, we see that in both cases the allowed regions for MSW solutions remain compara- 
ble in size to what they are calculated to be with realistic estimates of the nuclear physics 
uncertainties. The reason that the allowed regions are not particularly sensitive to the un- 
certainties in any one nuclear parameter is that there are a number of roughly comparable 
uncertainties from different input parameters pi . 



We conclude that one cannot greatly increase the accuracy with which MSW parameters 
can be determined by greatly decreasing the uncertainty in any single input parameter. 
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3. Sterile neutrinos 

How are the results given above changed if the oscillations involve sterile neutrinos? 

Figure [| shows that if the oscillations are between an electron-type neutrino and a sterile 
neutrino then the LMA and the LOW solutions are not allowed at the 99% C.L. ; only the 
SMA solution is possible . The LMA and LOW solutions are ruled out at a very high 
C.L.; the value of x 2 (min) = 19.0 (17.0) for the LMA (LOW) solutions. However, the SMA 
solution is still allowed, with x 2 (min) = 1.7 (which is the same as for the best-case with 
active neutrinos) and the best-fit solution is 

Am 2 = 4.3 x 10" 6 eV 2 , (8a) 



sin 2 20 = 6.9 x 10" 3 . (8b) 



The mass difference for the solution involving sterile neutrinos, Eq. (pa]) , is slightly smaller 
than the mass difference found for oscillations between active neutrinos (cf. Eq. Hal), but 
the mixing angle for sterile neutrinos, Eq. is slightly larger than for active neutrinos (cf. 
Eq.gb]). 

The suppression of the neutrino fluxes in the case of SMA conversion into sterile neutrinos 
is somewhat similar to the arbitrary reductions of the neutrino fluxes considered in Sec. |11 B 



The flux of pp neutrinos is unsuppressed, the 7 Be neutrino flux as well as other fluxes 
of intermediate energies are strongly suppressed, and the 8 B neutrino flux is moderately 
suppressed. The fit to the solar neutrino experiments using oscillations into sterile neutrinos 
is much better than the fit to the experiments made using an arbitrary linear combination of 
undistorted solar neutrino fluxes (cf. Sec. |11B| ). The main difference between the suppression 



due to sterile neutrinos and the suppressions considered in Sec. [II B| is that oscillations 
into sterile neutrinos suppress the lower energy part of the 8 B neutrino energy spectrum. 
This result is an illustration of our claim made in Sec. 11 CI that the evidence from the 



total experimental rates suggests the existence of a distortion of the solar neutrino energy 
spectrum. 

What should one expect for BOREXINO if oscillations involve sterile neutrinos? The 
allowed range for the 7 Be electron scattering rate is 

(0 cr ) 7 Be sterile _ n nnn +0.244 



BP98 



0.009^- (9) 



If oscillations occur to sterile neutrinos, the rate observed in BOREXINO will be — for almost 
the entire range of allowed parameters — less than the lowest expected rate if oscillations 
occur to active neutrinos (cf. Eq. |8| ). 

If the Caltech normalization for the 8 B production cross section (Sn = 22.4 kev b) is 
used instead of the INT normalization, the allowed regions are shifted only slightly from 
what they are in Figure |2] and Figure |j. The best estimate solutions for 5m 2 are changed 
by less than 7% relative to the values given in Eqs. flU)-® and Eq. (H). 
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B. Vacuum neutrino oscillations 



Figure |5| shows for vacuum neutrino oscillations the broad region of solutions allowed at 
99% C.L. . The best-fit vacuum solution is 

Am 2 = 8.0 x 10- n eV 2 , (10a) 
sin 2 29 = 0.75, (10b) 

which has a xLin = 4.3. 

What will BOREXINO observe if vacuum neutrino oscillations occur? The allowed 
range for vacuum oscillations of the 7 Be line is very large if the only constraints imposed are 
consistent with the total observed rates. We find at 99% C.L. that 

(0°") 7 Be vacuum _ n , fi +0.46 



BP98 



0.46!£S • (11) 



There is no allowed solution at the 99.7% C.L. for vacuum neutrino oscillations between 
u e and a sterile neutrino. The solution involving sterile neutrinos has % 2 (min) = 12.0. 

C. Arbitrary 8 B neutrino flux 



The value of the 8 B neutrino flux calculated in the standard solar model (cf. Ref. |20|| ) is 



more uncertain (+19% and —14%, la) than any of the other experimentally-important solar 
neutrino fluxes. It is therefore useful to consider what constraints are placed upon neutrino 



physics if the B flux is treated as a free parameter [ 48,33] . 



1. MSW solutions 

Figure |B| shows the allowed parameter space for MSW oscillations when the 8 B flux is 
allowed to take on arbitrary values. The best-fit SMA solution is 

Am 2 = 5.0 x 10~ 6 eV 2 , (12a) 
sin 2 20 = 3.5 x 10" 3 , (12b) 

which has Xmm = 0.86. The best-fit for the LMA solution is 

Am 2 = 1.6 x 10~ 5 eV 2 , (13a) 
sin 2 29 = 0.57, (13b) 

which has Xmin = 0.91. The LOW solution, 

Am 2 = 7.9 x 10" 8 eV 2 , (14a) 
sin 2 29 = 0.95, (14b) 
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has a much larger xLm = 7.2. For all three classes of MSW solutions the squared mass 
differences are changed by less than 20% if the 8 B neutrino flux is treated as a free parameter 
(cf. Eq. |j- Eq. |6|), although the values of sin 2 29 are changed by much large factors (a factor 
of two for the SMA solution) . 

Do the solar neutrino observations place useful limits on the value of the cross section 
factor, 5*i7, for the production of 8 B? 

Figure [7| shows the minimum value of y 2 obtained using as constraints the rates of the 
solar neutrino experiments but allowing the value of the 8 B flux to take on arbitrary values. 
The allowed range at 99% C.L. is: 0.4 < 5 , i 7 /5 , (INT) 17 < 2.0, which is much broader than 
the range allowed |^TJ by direct laboratory experiments. 

If oscillations between sterile neutrinos are considered, the constraints set on Sn by 
varying the solar 8 B flux and fitting to the solar neutrino data are again much larger than 
the uncertainty in the laboratory measurement of Syj (cf. Ref ||49||). 



2. Vacuum oscillations 

How are the inferred parameters for vacuum neutrino oscillations affected by permitting 
arbitrary values for the assumed 8 B neutrino flux, cf. Ref. [PJ? 

Figure § shows the expanded solution space that is allowed for vacuum oscillations if 
the 8 B neutrino flux is unconstrained. Instead of the rather limited parameter space that is 
permitted with the standard solar model value of the 8 B flux (see Fig [J), the squared neutrino 
mass difference can span the entire two order of magnitude range from Am 2 = 4 x 10~ 12 eV 2 
to Am 2 = 5 x 10~ 10 eV 2 if Sn is allowed to vary arbitrarily. The allowed range of the 8 B 
cross section factor at 99% C.L. is, for vacuum oscillations, 0.4 < Si 7 /S(INT)i 7 < 2.0. 

For completeness, we note that the best-fit values for the vacuum solution and arbitrary 
8 B flux are 

Am 2 = 8.4 x 10" n eV 2 , (15a) 



sin 2 2# = 0.98, (15b) 

which has a shallow Xmm = 0.94. The minimum occurs for a 8 B flux that is 1.9 times the 
value of the predicted BP98 8 B flux. 

This allowed range of Si 7 is very similar for both the MSW and the vacuum solutions. In 
both cases, the measurements by SuperKamiokande determine the upper and lower limits 
on the allowed 8 B flux. The lower limit corresponds to very little conversion from v e to 
other neutrino types. The lowest allowed value is obtained by reducing the best-estimate 
SuperKamiokande flux by the 3cr experimental uncertainty. Since the experimental uncer- 
tainty is small (lcr exp = 3%), the lower limit is close to the SuperKamiokande best-estimate 
value. The upper limit is much larger because essentially all of the 8 B flux observed by 
SuperKamiokande is or v T . Since the cross section [18| for electron- neutrino scattering 



by or v T is about six times smaller than the scattering cross section for z/ e , and and v T 
are not detected in the chlorine and gallium experiments, the upper limit for the total flux 
can be much larger (a factor of five in practice) than the observed flux. 
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D. Energy-independent suppression of neutrino fluxes 



Many authors have considered particle-physics models which predict energy-independent 
reductions of the solar neutrino fluxes (see, e. g., Ref. [|50|-|54|1 and references cited therein). 



In these scenarios, all neutrino fluxes are reduced by some particle-physics mechanism by 
exactly the same factor. One can test the goodness-of-fit of such scenarios by calculating the 



minimum x 2 for this case [55|. Using the predictions and the uncertainties in the BP98 solar 
model and the observed rates and their uncertainties given in Table [I], we have calculated 
the minimum % 2 for different neutrino oscillation scenarios. 

For arbitrary reduction factors a, the best-fit value is a = 0.48 for oscillation into active 
neutrinos and a = 0.50 for sterile neutrinos. The minimum Xmm = 12.0 for active neutrinos 
and Xmin = 19-3 for sterile neutrinos. Since there are two degrees of freedom in these 
cases, energy-independent oscillation into active neutrinos is ruled out at the 99.8% C.L. 
and energy-independent oscillation into sterile neutrinos is ruled out at the 99.99%C.L. . 
These results are changed only slightly (strengthened slightly) if the 8 B flux is allowed to 
vary independently of all other fluxes. In this case, we obtain for active neutrinos (sterile 
neutrinos) Xmm = H-0 (Xmin = 19.3) for 1 d.o.f., which is excluded at the 99.8% C.L. 
(99.999% C.L.) . 



For the model considered in Ref. |5T|, a = 5/9. In this case, Xmm = 14.3 (3 d.o.f.), and 



the model is ruled out at the 99.8% C. L. 



E. Summary of analysis of rates 

The principal results of this section are displayed in Fig. 0-Fig. except for the BOREX- 
INO predictions which are given in Eq. (]?]), Eq. @, and Eq. ([TT|) , and the demonstration in 
Sec. II V Dlthat constant-reduction solutions are disfavored. For the reader's convenience, we 



also present the numerical values of the best-fit neutrino mixing angles and mass differences 
for each case considered. 

In all of the calculations shown in the figures, we have used the INT normalization for 
the 8 B production cross section, Srr- However, this cross section is relatively poorly known 
and, unfortunately, at present it is largely a matter of judgment as to which normalization, 
INT or CIT, is most appropriate. We have therefore presented results for both normaliza- 
tions; the difference between neutrino parameters determined with the INT and the CIT 
normalizations is an indication of the magnitude of the uncertainties in neutrino parameters 
that are caused by uncertainties in the basic nuclear physics data. 

Particle physics solutions in which the suppression of the electron neutrino flux is inde- 
pendent of the neutrino energy are ruled out at more than the 99% C. L. for oscillations 



into both active and sterile neutrinos (see Sec. |I V D| ) . 

Table |TJ presents the minimum x 2 for the different neutrino oscillation scenarios that we 
have considered in this section. The primary entries in the table represent x 2 (min) with the 
INT normalization for Sn; the entries enclosed in parentheses were computed using the CIT 
normalization for Sn. The INT and CIT normalizations give about equally good fits to the 
observed rates in solar neutrino experiments, with the INT normalization slightly favoring 
the SMA MSW solution and the CIT normalization slightly favoring the LMA solution (see 
also Fig. 0). Indeed, the larger 8 B flux, the CIT normalization, implies a stronger suppression 
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by oscillations at higher energies, and consequently, a flatter dependence of the suppression 
factor on energy than is implied by the INT normalization. 



V. ZENITH- ANGLE DEPENDENCE OF RATES 

If MSW oscillations occur, the observed event rates in solar neutrino experiments can 
depend upon which region, if any, of the earth that the neutrinos traverse before reaching 
the detector. At times when the sun is below the horizon, or v T coming from the sun can 
be re-converted in the earth into the more easily detected v e . Thus for a certain range of 
neutrino parameters the observed event rate will depend upon the zenith angle of the sun, 
i.e., the angle that the sun makes with respect to the direction of the zenith at the position 



of the neutrino detector. This is known as the earth regeneration effect [56-58 . 



A. Expectations for survival probabilities 

Figure |9] shows the computed survival probabilities for electron type neutrinos as a func- 
tion of energy for the day (no regeneration), the night (with regeneration), and the annual 
average. The survival probabilities computed here for the best-fit points differ relatively 
little from the values calculated earlier (see Fig. 4 of Ref. fl43"|), despite the fact that in the 



present study we have taken account of the measured SuperKamiokande rate, the somewhat 
different fluxes from the BP98 model, the improved estimates for nuclear fusion cross sec- 
tions, and other updated data. The changes in the best-fit oscillation parameters caused 
by the use of the new data are relatively small (see Sec. |TV| ). The principal difference is 
that with the present parameters the expected Day-Night difference for the SMA solution 
is extremely small, not even visible in Fig. []. Note, however, that a potentially detectable 
Day-Night difference is predicted for the LMA solution at the higher neutrino energies. The 
Day-Night difference is expected to be relatively large for the LOW solution, but only at 
energies below 1 MeV. The Day-Night difference predicted by the LOW solution would be 
detectable by the BOREXINO experiment, which will observe the 0.86 MeV neutrino line 
from 7 Be electron capture. 

B. SuperKamiokande result 

The SuperKamiokande Collaboration has given a preliminary value for the difference 
between the event rates at night and during the day. After only 504 days of data, they 
obtain an initial estimate for this Day- Night asymmetry, A, of 

, D-N . , 

A = D + N = -0.023 ± 0.020(stat) ± 0.014(syst). (16) 

This estimate applies for events in which the recoil electron has an energy of at least 6.5 
MeV. The difference shown in Eq. (|16|) is in the direction that would be expected from 
regeneration in the earth (the sun is apparently brighter at night in neutrinos), but is small 
and is not statistically significant. 
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The SuperKamiokande Collaboration has also given ||,|10] a preliminary distribution of 
the event rates versus zenith angle in which the rates and their errors are plotted in 10 
angular bins. 



C. Constraints on neutrino parameters 

We have determined the constraints placed upon neutrino mixing parameters using the 
preliminary Day-Night asymmetry (Eq. |l(J) and also the 10-bin zenith angle distribution. We 
use the techniques and results contained in our previous detailed discussion of the expected 
zenith-angle dependence of solar event rates to analyze the 10 bin measurements of the 
zenith angle distribution reported by the SuperKamiokande Collaboration for the first 504 
days of observations JTOj . 

Figure |lO| a shows the results of a combined MSW \ 2 fit with active neutrinos of the 
SuperKamiokande zenith angle dependence and the total rates of the Homestake, GALLEX, 
SAGE, and SuperKamiokande experiments. The region in the figure that is excluded by the 
zenith angle dependence is shaded lightly and, almost touches the darkly shaded LMA and 



SMA allowed regions for the combined zenith angle and rate constraints. Figure |10| b shows 
the combined solutions using the SuperKamiokande Day-Night asymmetry and the four 
total rates. As shown previously from simulated data E|, the Day- Night asymmetry is 



more sensitive to the LMA solution and the zenith-angle distribution is more sensitive to 
the SMA solution. We see from Fig. |T0| b that the observations by SuperKamiokande of the 
Day-Night asymmetry have eliminated a sizable fraction, almost a half, of the parameter 
space for the LMA solution that is allowed if one only considers the total rates (cf. Fig. ^|). 
A small part (with sin 2 29 > 10~ 2 ) of the parameter space of the SMA solutions is excluded 
by the measured zenith-angle dependence. 



Figure |TT| a and Figure pJ] b show the analogous regions for MSW oscillations into sterile 
neutrinos. The only MSW solution for oscillation into sterile neutrinos is the SMA. The 
regeneration effect is not as effective a test of the sterile neutrino oscillation solution as it is 



for oscillation into active neutrinos (cf. Fig. ITT] and Fig. 



VI. PRELIMINARY SPECTRAL SHAPE 

If minimal standard electroweak theory is correct, the shape of the 8 B neutrino energy 
spectrum is independent of all solar influences to an accuracy of 1 part in 10 5 J55 |. The 



shape of the neutrino spectrum determines the shape of the recoil electron energy spectrum 
produced by neutrino-electron scattering in the detector. Therefore, any departure of the 
observed electron recoil energy spectrum from the shape predicted by standard electroweak 
theory would be a "smoking gun" indication of new physics. In this section, we compare 
the preliminary recoil electron energy spectrum reported by SuperKamiokande |T(J with 



the results calculated using the standard (undistorted) 8 B |1| neutrino spectrum. We 
also compare the observed energy spectrum with recoil energy spectra calculated assuming 
MSW or vacuum neutrino oscillations and use these results to constrain the allowed neutrino 
parameter space. 
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A. SuperKamiokande energy spectrum 

The SuperKamiokande Collaboration has made available |Tt| preliminary data, including 
estimated statistical and systematic uncertainties, from 504 days of operation in which 
the recoil energy spectrum is divided into 16 bins, with 15 bins having a width of 0.5 
MeV starting at 6.5 MeV and continuing to 14 MeV. The final bin includes events with 
energies from 14 MeV to 20 MeV and contains more counts than would be expected from 
an undistorted 8 B spectrum normalized at lower energies. 



B. Fit to undistorted energy spectrum 

The minimum \ 2 f° r the fit of the undistorted energy spectrum to the measured energy 
spectrum is 31 for 15 D.O.F. . The fit with the undistorted spectrum is acceptable only at 
slightly less than the 1% C.L., consistent with the results reported by the SuperKamiokande 
Collaboration at Neutrino 98 fllQfl. 



C. Fits to MSW and vacuum oscillations 



The simplest test |K) for a deviation from the standard recoil energy spectrum is to 



investigate whether the ratio, R, of the observed to the standard energy spectrum is a 
constant, which would be expected in the absence of a distortion. We have therefore fit 
the ratio R to a linear function using the measured ||10[ number of events and their quoted 
uncertainties in the 16 bins (from 6.5 MeV to 20 MeV) and the standard electron recoil 
energy spectrum []nj modified by the energy resolution of the SuperKamiokande detector. 
Thus 

R = R + S * (W e - 10 MeV), (17) 

where R represents the average event rate and So is the average slope that measures the 
deviation of the recoil electron energy spectrum from the undistorted shape. The total 
electron energy is W e . The energy spectrum of recoil electrons is determined by convolving 
the neutrino spectrum [[H]] with the calculated survival probability, the neutrino-electron 
scattering cross-sections fli~8" |, and the energy resolution function ||10|| . These integrations 



smooth the effect of distortions; therefore, the expected distortion can be well described by 
a simple linear dependence. 

Figure [12] shows the result of this calculation. The la, 2a, and 3a allowed regions are 
shown in the figure. We have taken account of the bin correlations between the systematic 
errors in the fashion explained in the Appendix. The best-fit point lies at Rq = 0.474 and 
S = 0.0153, with xLin = 23.5 for 14 d.o.f.. Therefore, a straight line is not a particularly 
good fit to the data; it is acceptable at the 5.3% C.L. . The main reason that the fit 
is acceptable only at a modest C.L. is that the last three bins indicate a deviation from 
a smooth extrapolation at lower energies. This deviation could be real or it could be a 
statistical fluctuation. 



Figure [12| also shows the ranges of the predicted slopes (S ) for different successful 



descriptions of the solar neutrino total experimental rates, i.e., the SMA active and sterile, 
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LMA, LOW, and vacuum oscillation solutions described in Sec. [TV]. All of the oscillation 
solutions shown in Fig. [12| permit a wide range of Rq that includes the range measured 
by the SuperKamiokande Collaboration. For visual convenience, we have not included the 
horizontal error bars, i. e., the uncertainties in R , for the neutrino oscillation solutions 
shown in Fig. |12|.p| Therefore, there are acceptable oscillation solutions of all four types that 
describe both the total rates in the four solar neutrino experiments and the recoil electron 
energy spectrum measured by SuperKamiokande. 

How powerful a constraint is the spectrum shape in determining the allowed neutrino 
oscillation parameters? 

Figure [13] shows, for both active and sterile neutrinos, the parameter region that is al- 
lowed by considering only the spectral information. Only a small fraction of MSW parameter 
space is consistent with the spectral data. It is instructive to compare the regions allowed 
by the rates only, shown in Fig. |2| and Fig. |], with the regions shown in Fig. |TJ, which are 
allowed by the spectrum data. The best-fit solutions considering only the spectrum, marked 
by dark points in Fig. [13], do not lie within the allowed regions for the global solutions 
discussed in Sec. |V11| . 

Figure [TJ] shows the region that is allowed for vacuum oscillations by the information 
from the SuperKamiokande energy spectrum. The dark circle shows the best-fit point. Just 
as for MSW oscillations, only a small fraction of the parameter space is allowed for vacuum 
oscillations by the spectrum constraint. The complementarity of the analysis that uses only 
the rates and the analysis that uses only the spectral data can be seen by comparing Fig. |5| 
and Fig. [TJ. 



VII. GLOBAL SOLUTION: RATES, ZENITH-ANGLE DISTRIBUTION, AND 

ENERGY SPECTRUM 

We discuss in this section simultaneous fits to all the available data. We consider solutions 
that describe the total rates in the four experiments (chlorine, SuperKamiokande, GALLEX, 
and SAGE) plus either the energy spectrum measured by SuperKamiokande or both the 
energy spectrum and the Day-Night asymmetry (or zenith-angle dependence) measured by 
S uper Kamiokande . 

We begin by showing in Sec. |V11 A| that an arbitrary combination of undistorted solar 
neutrino energy spectra is, independent of any astrophysical considerations, ruled out at the 
3.5cx level. We next show in Sec. |V11B| that the best global fit for the MSW solutions is 
acceptable at the 7% C.L. for active neutrinos (8% C.L. for sterile neutrinos). The best-fit 
global solution is very close the best-fit SMA solution when only the total experimental rates 
are considered. Finally, we show in Sec. [VII C| that the best global fit for vacuum oscillations 
is acceptable at the 5% C.L. . 



*The allowed range, cf. Sec. |fVj of R is 0.38 to 0.86 for the SMA active solution, 0.33 to 0.84 for 
the SMA sterile neutrin solution, 0.31 to 0.55 for the LMA solution, and 0.42 to 0.53 for the LOW 
solution. 
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The results of the simultaneous fits to all of the available data of the neutrino predictions 
are shown in Fig. |l5|-Fig. 17, which present the allowed regions for the different oscillation 
scenarios. 



A. Global fits: arbitrary undistorted fluxes 



What is the best fit to the total rates plus the SuperKamiokande spectrum and Day- 
Night asymmetry if we allow arbitrary values, subject only to the luminosity constraint |pT| . 
for the pp, 7 Be, 8 B, and CNO fluxes ? The minimum y 2 satisfies 

Xminimum( 3 rates + spectrum + D/N; arbitrary pp, 7 Be, 8 B, 13 N, 15 0) = 39.2 , (18) 

for 15 d.o.f.| This result is excluded at the 99.94% C.L., which corresponds to a 3.5<r 
deviation from minimal electroweak theory. 



B. Global fits: MSW solutions 



Fig. |l5| a displays the results of imposing on the MSW solutions the combined constraints 
of the total rates and the SuperKamiokande spectrum. The best-fit MSW solution con- 
sidering both total rates in the four experiments and the electron recoil energy spectrum 
measured by SuperKamiokande has 

Am 2 = 5.4 x 10~ 6 eV 2 , (19a) 



sin 2 20 = 6.3 x 1(T 3 . (19b) 

The minimum xLm = 26.5, which is acceptable at the 7% C.L. (for 17 d.o.f.), not a very good 
fit. The best-fit solution is very close to the SMA solution when only rates are considered 
(cf. Eq. f|). Including the spectrum in addition to the total rates, eliminates (at 99% C.L.) 
the entire LOW solution and a large portion of the LMA solution at higher Am 2 . The 
inclusion of the spectral constraint also eliminates for the SMA solution part of the region 
at smaller mixing angles that is allowed if only the total rates are considered. 

Fig. ^5] b shows the allowed region for MSW parameters when the constraints from the 
SuperKamiokande zenith-angle distribution is included together with the constraints from 
the four measured total rates and the SuperKamiokande electron recoil energy spectrum.Q 



§There are 3 d.o.f associated with the rates, 15 d.o.f associated with the spectra shape which has 
one overall normalization parameter that is variable, and 1 d.o.f. for the Day-Night asymmetry. 
All 5 fluxes are allowed to vary freely, subject only to the luminosity constraint. 

**We have not shown the combined fit including the Day-Night asymmetry, the spectrum con- 
straint, and the total experimental rates since the excluded region is very similar to what appears 
in Fig. p~5|b. We chose to display the fit with the zenith-angle distribution since this provides the 
best restriction on the SMA solution. 
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The best-fit solution is almost identical to what is obtained for the rates only case and for 
the case of rates plus zenith- angle constraint, namely, 

Am 2 = 5.4 x l(T 6 eV 2 , (20a) 

sin 2 20 = 5.5 x KT 3 . (20b) 

The minimum xLin = 37.2, which again is acceptable at the 7% C.L. (for 26 d.o.f.). Only the 
SMA solution survives at the 99% C.L. when the zenith-angle and the spectrum constraints 
are added to the constraints of the total experimental rates. The LMA and the LOW solution 
are marginally ruled out (Xmin = 47(49) for the LMA (LOW) solution). 

The Day-Night asymmetry parameter, A, which is defined by Eq. ([16]), lies in the range 



- 0.0048 < A < 0.025 , (21) 

for active neutrinos when all three sets of constraints are applied. The best-fit solution 
predicts a Day-Night asymmetry whose absolute value is less than 0.1%. 



How is the global solution changed if oscillations involve sterile neutrinos? Figure 16 
shows the allowed solution space when the measured total rates, energy spectrum, and 
zenith-angle dependence are all included as constraints. The best-fit solution is very close 
to what is obtained from the rates only (or the rates plus spectrum constraints) for sterile 
neutrinos. We find 

Am 2 = 4.0 x 10~ 6 eV 2 , (22a) 

sin 2 20 = 6.9 x 10~ 3 . (22b) 

The minimum xLm = 36.5, which is acceptable at the 8% C.L. (for 26 d.o.f.), slightly better 
than for active neutrinos. For sterile neutrinos, the Day-Night asymmetry parameter, A, 
lies in the range 

- 0.0031 < A stcrilc < 0.007 . (23) 

The best-fit solution predicts a value for the asymmetry parameter whose absolute value is 
less than 0.1%. 



C. Global fits: vacuum oscillations 



Fig. [17] shows the allowed regions for the vacuum oscillation parameters when the con- 
straints from the rates, the spectrum shape, and the Day-Night asymmetry are all included. 
The best-fit vacuum oscillation solution has 

Am 2 = 6.5 x 10~ n eV 2 , (24a) 

sin 2 29 = 0.75 . (24b) 

The minimum Xmm = 28.4, which is acceptable at the 6% C.L. (for 18 d.o.f.). This value 
of Xmm i s only slightly below the value of \ 2 = 30.5 that is found for the best-fit point in 
the 'spectrum-only' analysis (cf. Fig. [14]). The inclusion of the spectral constraints reduces 
considerably the domain of allowed solutions for vacuum oscillations (cf. Fig. [17] and Fig. ||) . 
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D. Global fits: predicted 7 Be v — e scattering rates 

We have calculated the range of allowed 7 Be neutrino-electron scattering rates (0.86 MeV 
line) that is consistent with all of the available solar neutrino data: total event rates, Day- 
Night asymmetry, and spectrum shape. The results are very similar to what was obtained 
earlier in the discussion of the rates only, see Eq. (0) , Eq. @ , and Eq. (]IlD . We conclude that 
the neutrino oscillation predictions for the BOREXINO experiment are robust. We find, in 
particular, that the global range predicted by the MSW solution with active neutrinos is 

(<M 7 Be SMA _ n O o+0.24 f r\ 
{(pa 7BP98 

The corresponding range for MSW oscillations into sterile neutrinos is 



(0C") 7 Be sterile _ n nnfi +0.25 



(0a> 



BP98 



0.00618;*,. (26) 



For vacuum oscillations, 



^) 7 Be vacuum _ n , r+0.33 



(0CT> 



BP98 



0.45^ . (27) 



The principal change caused by the imposition of the additional constraints due to the 
measured spectrum shape and Day-Night asymmetry is a modest shrinking of the allowed 
range for the vacuum solutions. 



E. Global fits: predicted T Be v e fluxes 



Raghavan |61[] has recently proposed a flavor-specific (neutrino-absorption) experiment 
that measures the total flux of 0.86 MeV v e reaching the earth from 7 Be electron captures 
in the sun. The preliminary name for this experiment is LENS. 

We have calculated the range of allowed 7 B v e survival probabilities for comparison 
with the results of a future LENS experiment. We again considered neutrino parameters 
consistent with all the available solar neutrino data: total event rates, Day-Night asymmetry, 
and spectrum shape. We give below the predicted values for the survival probability, P, 
for an electron-type 0.86 MeV neutrino created in sun to remain an electron-type neutrino 
when it reaches the terrestrial target. For MSW oscillations into active neutrinos we find 

-P 7 Be sma = 0.02^q;q5- (28) 
The corresponding range for MSW oscillations into sterile neutrinos is 

-P 7 Be sterile = 0.006^q;oq 2 . (29) 

For vacuum oscillations, 

-P 7 Be vacuum = 0.30^q 14 . (30) 

With current knowledge, the allowed range of P varies all the way from 0.00 to 0.72, 
which emphasizes the importance of this proposed measurement. 
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F. Global fits: energy-independent suppression 



The preliminary results of the SuperKamiokande measurement of the electron recoil 
energy spectrum provide additional evidence, beyond that available from just the total 
rates discussed in Section |IV D| , regarding the hypothesis [p r Q|-|5^] for an energy-independent 
suppression of the solar neutrino fluxes (but see also the discussion of Figure in Sec- 
tion VIII E|) f^|. Assuming an energy-independent suppression of the neutrino fluxes in the 
BP98 standard solar model [20], x min f° r oscillations into active neutrinos is 43 (17 d.o.f), i.e., 
this solution is disfavored at 99.95% CL. . If the 8 B neutrino flux is allowed to vary as a free 
parameter, x min = 42 ( 16 d.o.f), which is disfavored at 99.96% CL. . Energy-independent 
oscillations into sterile neutrinos are disfavored at 99.998% CL. (minimum x 2 — 50.3) if 
one allows an arbitrary 8 B neutrino flux. Maximally-mixed oscillations into active neutrinos 
which predict [^TJ a constant survival probability of 5/9 give Xmm = 45 (18 d.o.f.), which 
is disfavored at 99.96 % CL. .For arbitrary 8 B neutrino flux, \ 2 = 44.5 (17 d.o.f), which is 
disfavored at 99.97% CL. . 



VIII. DISCUSSIONS AND CONCLUSIONS 

We summarize and discuss in this section our principal conclusions. In Sec. |V111 A| , we 
discuss the indications that solar neutrino experiments are suggesting new physics. We 
present in Sec. |V111B| the robust inferences based only on the total event rates and in 
Sec. [VIII C| the implications of the absence of a statistically significant zenith-angle depen- 
dence for the SuperKamiokande event rates. In Sec. [VLI [D|, we analyze the strong constraints 
imposed by the electron recoil energy spectrum measured by SuperKamiokande. Our global 



analysis is described in Sec. [VIII E| , where we present the implications of the combined con- 
straints from the total rates, the zenith-angle dependence, and the spectrum shape. We 
summarize in Sec. [VIII F| our overall view of the solar neutrino situation. 

For a given hypothesis, MSW or vacuum oscillations (active or sterile neutrinos), we 
search numerically for the best-fit solution and quote the CL. for acceptance or exclusion 
based upon the relevant x min and the appropriate number of d.o.f. . The acceptance levels 
found in this paper seem somewhat higher than reported by the SuperKamiokande Collab- 
oration in their review talk at Neutrino 98, but we think that this difference is largely due 
to the fact that we have searched for the best-fit solution rather than test the acceptability 
of previously recognized solutions and that we have used the BP98 rather than the BP95 
predictions in analyzing the total rates. To the best of our knowledge, the present paper is 
the first to determine global solutions that take account of the measured total rates in all of 
the solar neutrino experiments, as well as the SuperKamiokande zenith-angle distribution 
and the SuperKamiokande electron recoil energy spectrum. 



^Different physical mechanisms that might lead to an energy-independent supression may predict 
different zenith angle dependences for the detector event rates. Therefore, we have not included 
information about the zenith-angle distribution in the global fits for energy-independent supression. 
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As described below in Sect. |VIIIE|, the best-fit global MSW solution for active neutrinos 



is: Am 2 
l(T 6 eV 2 
6.5 x 10" 



, sin 

-li 



5.4 x 10 6 eV 2 , sin 2 2# = 5.5 x 10 3 (and for sterile neutrinos: Am 2 
2 26 



eV 2 



, sin 



6.9 x 10~ 3 
2 26 = 0.75 



For vacuum oscillations, the best-fit solution is: 



= 4.0 x 
Am 2 = 



A. Indications of new physics 

The results from the first phase of the SuperKamiokande experiment have strengthened 
the inference that new physics is required to describe solar neutrino experiments. 

If we consider only the total rates in the solar neutrino experiments, then there is no 
linear combination of the undistorted neutrino energy spectra that can fit the available data 
at the 3cr level (see the discussion in Sec. [11 B 1| and Sec. [11 B 2| ). This result, whose physical 



basis is described in Sec. [II Q , is independent of any astrophysical arguments regarding the 
basic correctness of the solar model. In particular, the 3a discrepancy ignores the additional 
evidence provided by helioseismological measurements, which agree to high precision (0.1% 
r.m.s in sound speeds, see Ref. [f^J) with the predictions of the standard solar model. The 
Standard Solar Model predictions are inconsistent with the observed rates in solar neutrino 



experiments at approximately the 20<r level (see Sec. 11 A| ). 



The data from the total rates alone indicate that the v e energy spectrum from the sun 
is distorted, i.e., the survival probability for electron type neutrinos to reach the earth is 
energy dependent. 

If we impose the additional constraints from the measured SuperKamiokande spectrum 
and Day-Night asymmetry as well as the total rates, then an arbitrary linear combination 
of the pp, 7 Be, 8 B, and CNO neutrino fluxes is ruled out at the 3.5cr confidence level (see 
Sec. |V11 A|). This minimum discrepancy again ignores all information about the solar model. 



B. Average event rates 

The most robust results of the solar neutrino experiments so far are the total observed 
rates, which are summarized in Table |. We have therefore evaluated accurately the al- 
lowed regions of neutrino parameters for either active or sterile MSW or vacuum neutrino 
oscillations, using the total rates as the only constraints. If neutrino oscillations are in- 
deed occurring, a subset of the parameters that are consistent with the total rates must 
also be consistent with the other measured quantities (electron recoil energy spectrum and 
Day-Night asymmetry), a proposition that we have also tested. 

The principal results considering only the total rates are displayed in Fig. ||-Fig. §. We 
have also given in Sec. ^ the best-fit neutrino oscillation parameters and mass differences 
for each scenario that we have discussed. 

The most important change in the allowed range of neutrino parameters compared to our 
previous work |43| (which was prior to the announcement of the SuperKamiokande results) 



is that the SMA solution is shifted to somewhat smaller mixing angles, from the earlier value 
of sin 2 26 = 8.7 x 10" 3 @ to the current-fit value of sin 2 26 = 6.0 x 10" 3 (Eq. ga|). The main 
causes of this shift are the smaller predicted 8 B neutrino flux (for the INT normalization) 
and the the lower SuperKamiokande rate (somewhat lower than the Kamiokande rate). If 
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the low energy cross section factor for the 7 Be(p, 7) 8 B reaction is treated as a free parameter, 
then the best-fit mixing angle for the SMA solution decreases to an even smaller value of 
sin 2 29 = 3.5 x 10~ 3 (Eq. |12a| ). Further improvements in the most important nuclear cross 
sections are unlikely to reduce significantly the size of the allowed neutrino parameter regions 
(cf. Fig. |]), although further moderate shifts in the best-fit values may be anticipated. 

The relatively modest shift in mixing angles from sin 2 29 = 8.7 x 10~ 3 to sin 2 29 = 6.0 x 
10 -3 has the effect of reducing drastically the predicted Day- Night asymmetry from what was 



expected earlier ||43|| , i.e., 1.8% for the best-fit SMA solution, to a probably un-measurably 
small 0.35% asymmetry for the current best-fit SMA solution (and is even smaller for the 
best-fit global solution, see Sect. |V111 F| below). The Day-Night asymmetry is obviously 
sensitive to details of the neutrino parameter solutions. 

Some theoretical models are already strongly disfavored by the constraints of the total 
rates alone. For example, models in which the suppression of electron type neutrinos is 
independent of energy are excluded at the 99.8% C. L. or more, depending upon the precise 
scenario considered (see the results described in Sec. |VD|). 

C. Zenith- Angle dependence of rates 

No statistically significant Day-Night asymmetry or zenith-angle dependence of the solar 
neutrino event rate has been detected so far by the SuperKamiokande experiment |8HT0 



The small-value of the observed Day-Night asymmetry excludes a large part of the LMA 
region in neutrino parameter space that is allowed if only the total solar neutrino rates are 
considered (see Fig. P~C|) . 

D. Electron recoil energy spectrum 

The electron recoil energy distribution reported by SuperKamiokande is inconsistent 
with no distortion at about the 99% C.L. . On the other hand, all of the popular neutrino 
oscillation solutions determined from the total rates only (SMA, LMA, and LOW MSW 
solutions and vacuum oscillations) provide acceptable, although not excellent, fits to the 
recoil energy spectrum (cf. Fig. [12]). The observed distortion of the spectrum, i.e., the overall 
slope parameter, is in the direction predicted by the SMA solution. However, as emphasized 
by the SuperKamiokande Collaboration |1(J, the fits using the oscillation solutions for active 



neutrinos that are preferred by the total rates are not particularly good. The last three 
points in the recoil spectrum are somewhat higher than would have been expected from an 
undistorted 8 B neutrino spectrum. 



Figure [13] shows the regions of MSW parameter space that are allowed for both active 
and sterile neutrinos by imposing only the constraint of consistency with the electron recoil 
energy distribution. The best-fit solutions for the spectrum only, marked by dark points in 
Fig. [H|, do not fall within the allowed regions determined by the global fit to all of the data 
(see Sec. |V111E| below). For vacuum oscillations, Fig. [H] shows the regions allowed by the 
spectral data alone. 
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E. Global analysis 



The combined constraints from the total rates, the zenith-angle dependence (or the Day- 
Night Asymmetry), and the electron recoil energy spectrum provide the most comprehensive 
test of neutrino oscillation descriptions of solar neutrino experiments. 

The allowed solution space for MSW oscillations, including the rates and the energy 
spectrum, is shown in Fig. 15 a and the allowed solution space including all three sets of 
constraints, rates, energy spectrum, and zenith-angle dependence, is shown in Fig. pT5|b. The 
best-fit solution, in both cases, is close to the best-fit SMA solution when only the total rates 
are considered. The best-fit MSW solutions shown in both Fig. |T5| a and Fig. |T5] p are not 
particularly good, but are acceptable at the 7% C. L. . When all three sets of constraints are 
imposed, the SMA solution space is reduced somewhat in size compared to what is allowed if 
only the total rates are considered. For sterile neutrinos, the SMA solution is acceptable at 
the 8% CL. ; the allowed region, which is similar to what is obtained for active neutrinos, is 



shown in Fig. [16] . The LMA and LOW solutions are excluded for active and sterile neutrinos 
at the 99% C. L. when all three constraints, rates, spectrum, and Day-Night asymmetry, 
are included. 

The global fits of MSW solutions to all the available data predict a Day-Night asymmetry, 
defined by Eq. (jTB]), for the SuperKamiokande experiment (total energy > 6.5 MeV) that 
lies in the range -0.0048 < A < 0.025 for active neutrinos and -0.0031 < A < 0.007 for 
sterile neutrinos. The best-fit solutions predict an asymmetry whose absolute value is less 
than 0.1% , un-measurably small, for both active and sterile neutrinos. 

For vacuum oscillations, Fig. [17| shows the allowed solution space when the rates, energy 
spectrum, and zenith-angle dependence are all imposed. The best-fit vacuum solution is 
acceptable at the 6% C. L. . The inclusion of the Day-Night and the spectral constraints 
reduces considerably the parameter regions for vacuum oscillations (cf. Fig. |5|). 

Figure [18] shows how well (or poorly) the calculated global oscillation solutions fit the 
observed SuperKamiokande |10[ electron recoil energy spectrum. Each of the three panels 
shows for a different global oscillation solution (described in Section |VI1|) the ratio of the 
number of electrons in a given energy bin to the number expected using the undistorted neu- 
trino energy spectrum [O] and the electroweak neutrino-electron scattering cross sections 



with radiative corrections JTB|. The data shown were presented by the SuperKamiokande 



collaboration at Neutrino98 |T0| . In computing the predictions for the different global oscil- 



lation solutions, we included the reported SuperKamiokande energy resolution and trigger 
efficiency function. All three of the best-fit global solutions fall well below the measured Ra- 
tio in the three highest energy bins. The possibility that this discrepancy at high energies 
might be due to a larger-than-expected cross section for the production of hep neutrinos 



is discussed extensively by Bahcall and Krastev f62|| . The global best-fit vacuum oscilla- 
tion solution shows a slight upturn at the lowest energies that might be detectable when 
SuperKamiokande obtains data with a threshold that extends to 5 MeV . 

The neutrino oscillation predictions for the neutrino-electron scattering rate of the 0.86 
MeV 7 Be line, which will be measured by the BOREXINO neutrino experiment, are robust. 
Very similar predictions are obtained using only the total rates as a constraint and using 
the presently available data on the spectral shape and Day-Night asymmetry as well as the 



total rates. The results are summarized in Eq. ([^), Eq. (p6|), and Eq. (|27|) 
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We have also evaluated the globally allowed range of the v e survival probability which 
could be compared with future measurements of the v e 7 Be neutrino flux, in a charged 
current (LENS) experiment. The results are given in Eq. (|28|), Eq. (p9|), and Eq. (|30|). 



F. Summation 

Different neutrino oscillation scenarios, including SM A MSW conversion to either active 
or sterile neutrinos and vacuum conversion to active neutrinos, give similar quality descrip- 
tions of the solar neutrino data, acceptable at the level 6%-8%. Moreover, with the strong 
evidence for atmospheric neutrino oscillations than is now available, an oscillation solution 
of the solar neutrino problem seems even more probable than before. In addition, any de- 
scription of solar neutrinos that does not include some new physics that causes an energy 
dependence in the v e survival probability is strongly disfavored by the combined data from 
solar neutrino experiments, while standard solar model calculations accurately predict the 
helioseismologically measured sound velocities in the sun. 

However, with the existing data it is not possible to determine which kind of transi- 
tion^) solar neutrinos undergo. Fits of neutrino oscillation scenarios to just the total event 
rates in the chlorine, Kamiokande, GALLEX, SAGE, and SuperKamiokande experiments, 
which may be the most robust currently available experimental information, suggest neu- 
trino parameters in which the Day-Night asymmetry and the spectral distortion are difficult 
to measure. The situation may improve significantly when data from the lower-energy spec- 
tral bins of SuperKamiokande, at 5.5 MeV and 6.0 MeV, are available. Crucial information 
will also be provided by the future SNO [0, BOREXINO JTJ, and GNO @ experiments. 



We suspect that the unique description of solar neutrino phenomena will require global 
analysis of all of the available data. We hope that this paper is a useful step along the path 
toward a complete solution of the solar neutrino problems. 
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APPENDIX: ANALYSIS DETAILS 

We describe here the methods used in our statistical analysis of the solar neutrino data. 



24 



1. Rates only 



In the analysis which includes only the event rates in the different detectors, we use the 
procedure described in detail in Ref. H34]] . The y 2 for the combined fit is defined as: 

X 2 (Rates) = S M=M (i?f - CTy^f - ( Al ) 

Here Rf 1 is the theoretically predicted event rate in the i-th detector (chlorine, SAGE, 
GALLEX and SuperKamiokande), which takes into account the electron neutrino survival 
probabilities, i?^ xp is the corresponding experimentally measured event rate, and V^j is the 
error matrix which is a function of the theoretical uncertainties (nuclear cross-sections, age, 
luminosity and heavy element abundances in the standard solar model) as well as the experi- 
mental statistical and systematic errors from each experiment (see Table |). The theoretical 



uncertainties are determined in BP98 [20|. Following BP98, we have included theoreti- 
cal uncertainties due to the diffusion of elements in the sun and to the electron capture 
cross section on 7 Be, not previously included as uncertainties in neutrino oscillation studies. 
A detailed description of how the individual uncertainties are calculated can be found at 
http: / /www.sns.ias.edu/^jnb/SNdata| in the menu item 'Solar neutrino rates, fluxes, and 



uncertainties.' Important correlations exist between the neutrino fluxes that must be in- 
cluded correctly in the evaluation of the error matrix or the calculated allowed regions will 
be incorrectly reduced in size. 



2. Recoil electron spectrum 

Next we describe the analysis of the recoil electron spectrum in SuperKamiokande which 
includes a computation of the following \ 2 function: 

X 2 (Spectrum) = S iJ=ljl6 (aS* h - Sf^Wy^aSf - Sj Cxp ). (A2) 

Here Sf 1 is the theoretically predicted event rate for the i-th energy bin in SuperKamiokande 
and S exp is the corresponding experimentally measured event rate. In this paper, we use 
the measured electron recoil energy spectrum data presented by the SuperKamiokande Col- 
laboration [|H| at Neutrino 98 and the recently determined undistorted spectrum jOJ. The 



statistical as well as the systematic experimental errors are included. Since the latter are 
fully correlated, we assume a correlation coefficient of 1 between each pair of bins. Neglect 
of these correlations would lead to an appreciable increase of the regions that is excluded. 
The entries in the covariance matrix are: 

TJ/ stat state , syst syst /AQ\ 

Wij — a { a j dij + a/ a/ . (AS) 



The coefficient a in Eq. (|A2|) is an overall normalization coefficient and is varied as a 
free parameter independent of Am 2 and sin 2 29. This variation reflects the fact that we are 
interested here in a test of the shape of the measured spectrum and not in the overall event 
rate in SuperKamiokande. 

We calculate the theoretically expected event rates in the individual bins including the 
quoted energy resolution and trigger efficiency function in SuperKamiokande. We neglect 
the uncertainty in the energy normalization since, according to the SuperKamiokande Col- 
laboration, it is less than 1%. 
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3. Day-Night asymmetry and zenith angular dependence 



The Day-Night asymmetry measured by SuperKamiokande [TO is not significantly dif- 
ferent from zero, which limits the allowed neutrino oscillation parameters. Assuming 3cr 
errors and combining the statistical and systematic errors in quadrature, the asymmetry 
A exp , defined by Eq. (jl6[), is 

- 0.096 < A exp < 0.050. (A4) 



The lower limit in Eq. (|A4|) is not especially useful since we have found by explicit 
calculation that asymmetries predicted by MSW parameters in the entire range 10" 9 eV 2 < 
Am 2 < 10 _3 eV 2 and 10 -4 < sin 2 29 < 1.0, never take on very large negative values (A > 
—0.005). However, the upper limit is exceeded for certain choices of previously-allowed 
oscillation parameters; these parameters are therefore excluded by the SuperKamiokande 
measurement. 

Since the significant theoretical uncertainties cancel in the ratio, we define \ 2 for the 
Day-Night asymmetry as: 

X 2 (D/N) = (A th - A cxp )/AA cxp , (A5) 

where A t ^ is the theoretically predicted asymmetry due to neutrino regeneration in the earth 
(zero for vacuum oscillations) and Ay4 exp = 0.024 is the combined statistical plus systematic 
error given by the SuperKamiokande Collaboration. 

Figure [TI] and Figure [TT] show the excluded regions for active and sterile neutrinos, 
respectively. Note that the relatively large uncertainties in the total 8 B neutrino flux cancel 
in forming the asymmetry ratio and therefore do not affect the excluded regions. Also, 
we have verified that at the present level of accuracy the excluded regions are not affected 
significantly by the uncertainty in the undistorted 8 B neutrino spectrum. 

The measured zenith angular distribution provides another constraint on neutrino oscil- 
lation parameters. We have calculated the predicted zenith angular distribution of events for 
many points in the sin 2 29 - Am 2 plane and have compared the predicted distributions with 
the binned zenith-angle distribution presented ||10| by the SuperKamiokande Collaboration 
at the Neutrino 98 conference. 

The x 2 m this case is defined by 

X 2 (Zenith) = S ij=1 , 10 (aZf - Zf^U^aZf - Zf p ). (A6) 

The predicted distribution for the zenith angle dependence was calculated by assuming 
uninterrupted operation of the SuperKamiokande detector during 504 days starting from 
April 1, 1996. This assumption is not precisely correct since the detector has been occa- 
sionally shut down for calibration and maintenance. However, since these interruptions are 
a small fraction of the data taking period, our calculation is not seriously degraded by the 
lack of published information about the detailed operating schedule. 

The SuperKamiokande Collaboration can perform a unique public service by testing the 
sensitivity of their inferences based upon the zenith angle dependence to the assumptions 
regarding the operating schedule. They can test at what quantitative level the small effects 
due to the precise operating schedule affect the conclusions regarding neutrino parameters. 
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The matrix U that appears in Eq. ( |A6|) is assumed, in our calculations, to be diagonal. 
Correlations between the 10 bins of the zenith angular distribution might conceivably arise 
from systematic effects (for example, electrons falling in different bins being detected by the 
same photomultiplier tubes). Although these correlations are expected to be small, it would 
be very useful if the SuperKamiokande Collaboration were to publish their estimates of any 
correlations that cause U to be non-diagonal. 



4. Straight line fit 

As discussed in Sec. [VI Q , the ratio, R, of the measured to the expected spectral energy 
distribution can be fit by a straight line, as in Eq. flj"7l . For an undistorted spectrum, R is a 
constant independent of energy. 

The x 2 i n this case is defined as a sum over the 16 energy bins of the electron spectrum: 

X 2 (Line) = S iij=1 , 16 (aL* h - SF)Wrf(aL* - Sf p ) , (A7) 

where Lf 1 = R + S (E eji — E ), E = 10 MeV is a conveniently-chosen energy, E e>i is the 
total electron energy in the i-th bin (we use the energy in the middle of each energy bin), 
and Sq and Rq are arbitrary parameters which are varied until the minimum x 2 is found. 
The covariance matrix W is the same as the covariance matrix defined in Eq. (|A2|). 



5. Global fits 

In Sec. [VI I| , we combine in different ways the constraints from the measured total rates, 
the spectrum shape, and the Day-Night asymmetry (or zenith-angle dependence). In each 
case, we add the x 2 f° r eacn data set ( rates, energy spectrum, and Day-Night asymmetry 
or zenith-angle distribution). The individual x 2 are defined in previous sections of this 
appendix and are treated as independent in the global fits. Thus in Sec. |V11| the effective 
X 2 are always the sum of two or three independent terms. 
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TABLES 



TABLE I. Solar neutrino data used in the analysis. The experimental results are given in SNU 
for all of the experiments except Kamiokande (and SuperKamiokande), for which the result is 
expressed as the measured 8 B flux above 7.5 MeV (6.5 MeV) in units of 10 cm - 2 s — 1 at the earth. 
The ratios of the measured values to the corresponding predictions in the Bahcall-Pinsonneault 



standard solar model (BP98) of Ref. [20| are also given. The INT normalization (or, in parentheses, 
the CIT normalization) is used in these calculations. Only experimental errors are included in 
the column labeled Result /Theory. The results cited for the Kamiokande and SuperKamiokande 
experiments assume that the shape of the 8 B neutrino spectrum is not affected by physics beyond 
the standard electroweak model. 

Experiment Result Theory Result/Theory Reference 



Homestake 


2.56 


±0.16 ±0.16 


7.7± 


1.2 
1.0 


(8.81JS) 


0.33 ±0.029 


(0.29) 


|451 


Kamiokande 


2.80 


±0.19 ±0.33 


5.15 


+1.0 
-0.7 


(6-ltoi) 


0.54 ±0.07 


(0.46) 


i 


GALLEX 


77 


5 ± 6.2 


129 


+8 
-6 


(i3ii 9 7 ) 


0.60 ± 0.06 


(0.59) 


I 


SAGE 




66.6 + H 


129 


+8 

-6 


(i3ii 9 7 ) 


0.52 ±0.06 


(0.51) 


@ 


SuperKamiokande 


2.44 


± ±0.05 IT, 


5.15 


+1.0 
-0.7 




0.474 ± 0.020 


(0.39) 






TABLE II. Minimum x 2 f°r different neutrino oscillation solutions of the solar neutrino prob- 
lems The reference solar model corresponds to the INT normalization (or, in parentheses, the 
CIT normalization) of <Si7(0). Results are given for oscillations into either either active or sterile 
neutrinos. 



Solution 




MSW 




Vacuum Oscillations 




SMA 


LMA 


LOW 




active 


1.7 (2.8) 


4.3 (1.9) 


7.4 (8.2) 


4.3 (2.7) 


sterile 


1.7 (1.6) 


19.0 (17.0) 


17.0 (16.0) 


12.0 (11.2) 
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FIG. 1. Predictions of standard solar models since 1988. The figure shows the predictions of 19 
standard solar models in the plane defined by the 7 Be and 8 B neutrino fluxes. The abbreviations 
that are used in the figure to identify different solar models are defined in the bibliographical item, 



Ref. |B4|. We include all standard solar models with which we are familiar that were published in 
refereed journals in the decade 1988-1998. All of the fluxes are normalized to the predictions of the 



Bahcall-Pinsonneault 98 solar model, BP98 [20|. The rectangular error box defines the 3a error 
range of the BP98 fluxes. The best-fit 7 Be neutrino flux is negative. At the 99% C.L., there is no 
solution with all positive neutrino fluxes if the fluxes of CNO neutrinos are arbitrarily set equal to 
zero. There is no solution at the 99.9% C.L. if the CNO neutrinos are fixed at their standard solar 
model values. All of the standard model solutions lie far from the best-fit solution, even far from 
the 3cr contour. 
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FIG. 2. MSW solutions: rates only. The figure shows the allowed regions (99% C.L.) in Am 2 
— sin 2 29 parameter space for the MSW solution. The best fit points are indicated by dark circles. 
Only the total event rates in the chlorine, SuperKamiokande, GALLEX, and SAGE experiments are 
considered; the solar neutrino data are summarized in Table |. The 8 B neutrino flux corresponds to 
the "INT normalization" of Si7(0). The neutrino transitions in the sun are assumed to be between 
active neutrinos {y e — > z/„ or v e — > u T ). 
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FIG. 3. Effect of uncertainties in 534 and Sn on the allowed MSW solution space. The panels 
are the same as in Fig. ^ except that in the top panel the uncertainty in £34 is set equal to zero 
and in the lower panel the uncertainty in Sn is set equal to zero. There is no large reduction in 
the MSW allowed regions if either £34 or £17 is presumed to be negligible. 
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FIG. 4. Sterile neutrinos: rates only. The figure shows the regions in Am 2 — sin 2 29 parameter 
space that the total rates in the chlorine, SuperKamiokande, GALLEX, and SAGE experiments 
allow for MSW oscillations between an electron type neutrino and a sterile neutrino (y e — > u s ). 
Other conditions are the same as for Figure 0. 
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FIG. 5. Vacuum oscillations: rates only. The figure shows the regions in Am 2 — sin 2 28 
parameter space that the total rates in the chlorine, SuperKamiokande, GALLEX, and SAGE 
experiments allow for vacuum neutrino oscillations between active neutrinos. The best-fit point is 
indicated by a dark circle. The experimental rates are summarized in Table |. 
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FIG. 6. Variable 8 B flux and MSW solutions. The figure shows the allowed regions at 99% 
C.L. in the Am 2 - sin 2 29 plane for the MSW solutions with an arbitrary 8 B neutrino flux (treated 
as a free parameter). 
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FIG. 7. Minimum x 2 (MSW solution) as a function of the boron neutrino flux. The oscillations 
are between active neutrinos [y e — > or v e — > f r ). The reference boron flux in the BP98 solar 
model corresponds to INT normalization of Sn(0). 
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FIG. 8. Allowed regions in Am 2 - sin 2 29 parameter space for vacuum oscillations with an 
arbitrary 8 B neutrino flux. The oscillations are assumed to occur between active neutrinos. 
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Survival Probabilities 




FIG. 9. Survival probabilities for MSW solutions. The figure presents the yearly-averaged 
survival probabilities for an electron neutrino that is created in the sun to remain an electron 
neutrino upon arrival at the SuperKamiokande detector. The best-fit MSW solutions including 



regeneration in the earth are described in Sec. IV. The full line refers to the average survival 



probabilities computed taking into account regeneration in the earth and the dotted line refers to 
calculations for the daytime that do not include regeneration. The dashed line includes regeneration 
at night. There are only slight differences between the computed regeneration probabilities for the 
detectors located at the positions of Super-Kamiokande, SNO and the Gran Sasso Underground 
Laboratory (see Ref. p3f ). 
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FIG. 10. Angular distribution exclusion. The region that is excluded by the SuperKamiokande 
zenith-angle distribution or by the SuperKamiokande Day-Night measurement is shown as the light 
shaded area in the upper and lower panels, respectively. The darker shaded regions in the upper 
panel are the regions that are allowed by the total measured rates from the chlorine, GALLEX, 
SAGE, and SuperKamiokande experiments plus the zenith angle distribution measured by Su- 
perKamiokande. The darker shaded regions in the lower panel are allowed by the four measured 
rates and the SuperKamiokande Day-Night asymmetry. 
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FIG. 11. Angular distribution exclusion for sterile neutrinos. This figure is the same as Fig. 1C 
except that the present figure refers to oscillations into sterile neutrinos. 
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FIG. 12. Deviation from an undistorted energy spectrum. The la, 2a, and 3a allowed regions 
are shown in the figure. The ratio of the observed counting rate as a function of electron recoil 
energy [PJ to the expected undistorted energy spectrum [fyl was fit to a linear function of energy, 



with intercept Rq and slope Sq (see Eq. 17). The five oscillation solutions discussed in Sec. IV 



SMA active and sterile, LMA, LOW, and vacuum oscillations, all provide acceptable fits to the 
data, although the fits are not particularly good, see text in Sec. VI. 
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FIG. 13. Spectrum shape: allowed region for MSW oscillations. Figure 13a refers to active 
neutrinos and Figure [I~3|b refers to sterile neutrinos. Each panel shows the region in MSW solution 
space that is allowed by the SuperKamiokande [^] measurements of the recoil energy spectrum 
from the scattering of 8 B neutrinos by electrons. The best-fit solution when only the energy 
spectrum is considered is shown by dark points in Figure 13a and Figure Ol 



It is interesting to 

compare the regions allowed by the spectral energy distribution with the regions allowed by the 
total rates, cf. this figure with Fig. ^ and Fig. |I[ 
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FIG. 14. Spectrum shape: allowed region for vacuum oscillations. The figure shows the allowed 
region in the parameter space of vacuum oscillations that is permitted by the SuperKamiokande [1C] 
measurements of the recoil electron energy spectrum, cf. Fig. ||. The dark point shows the best-fit 
point considering the measured spectrum as the only constraint; the value of Xmm = 30.5. 
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FIG. 15. Global fits: MSW solutions. Figure |15| a shows the regions in MSW parameter space 
that are consistent with the total rates observed in the four solar neutrino experiments (chlorine, 
SuperKamiokande, GALLEX, and SAGE) and the measured SuperKamiokande electron recoil 



energy spectrum. Figure 15b shows the only allowed region in MSW parameter space that is 
consistent with the combined constraints from the four measured rates and the electron recoil 
energy spectrum and zenith angle distribution that are measured by SuperKamiokande. Contours 
are drawn at 99% C.L. 
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FIG. 16. Global fits: sterile neutrinos. The figure shows the allowed parameter region for MSW 
oscillations into sterile neutrinos that is consistent with the measured total rates, the zenith-angle 
distribution, and the recoil electron energy spectrum. Contours are drawn at 99% C.L. . 
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FIG. 17. Global fits: vacuum solutions. The figure shows the allowed parameter region for 
vacuum oscillations that is consistent with the measured total rates, the recoil electron energy 
spectrum, and the Day-Night asymmetry. Contours are drawn at 99% C.L. . 
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FIG. 18. Global best fits versus measured energy spectrum. The three panels compare the 



global neutrino oscillation solutions discussed in Section VII versus the electron energy spectrum 
measured by SuperKamiokande and reported at Neutrino98 [ 10 1 . The quantity, Ratio, that is 
plotted is the ratio of the number of electrons in a given energy bin, E, to the number that is 



calculated using the standard, undistorted B neutrino energy spectrum [19| and electroweak neu- 
trino-electron scattering cross sections with radiative corrections pl|. The no oscillation solution 



is a horizontal line, which, following the SuperKamiokande collaboration [10], is normalized to the 
BP95 prediction [27] and lies at Ratio = 0.37. 
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